Realization of 2D Spin-orbit Interaction and Exotic Topological Orders in Cold Atoms 
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Majorana zero bound mode exists in the vortex core of a chiral p+ip superconductor or superfluid, 
which can be driven from an s-wave pairing state by two-dimensional (2D) spin-orbit (SO) coupling. 
However, a 2D Rashba-type SO interaction is not experimentally realistic for an atom gas in the 
continuum. We propose here a novel scheme based on realistic cold atom platforms to generate 2D 
SO interactions in a blue-detuned square optical lattice, and predict both the quantum anomalous 
Hall effect and chiral topological superfluid phase in the experimentally accessible parameter regimes. 
This work opens a new direction with remarkable experimental feasibility to observe non-Abelian 
topological orders in cold atom systems. 
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Introduction.— The search for non-Abelian Majorana 
fermions has been a focus of both theoretical and exper- 
imental studies in condensed matter physics, driven by 
both the pursuit of exotic fundamental physics and the 
applications in fault-tolerant topological quantum com- 
putation [l-SJ. Majorana zero modes (MZMs) are pre- 
dicted to exist in the vortex core of a two-dimensional 
(2D) intrinsic (p -I- zp)-wave SC and at the ends of a ID 
p-wave SC [T]. Recent studies show that effective ID 
and 2D p-wave pairing states can be driven from s-wave 
SCs by spin-orbit (SO) interactions in the heterostruc- 
tures formed by conventional s-wave SCs and topolog- 
ical insulators [4 or semiconductors |5l [6], leading to 
MZMs in the case of odd number of subbands cross- 
ing the Fermi level. Following the theoretical prediction, 
Majorana end states have been suggestively observed 
in the semiconductor nanowire/ s-wave superconductor 
heterostructures [7] through tunneling transport mea- 
surements 8 . However, demonstrating the non-Abelian 
statistics of MZMs in braiding operation [2] , as an unam- 
biguous verification, is a far more demanding task and 
not yet available in solid-state experiments. 

On the other hand, the recent great advancement in re- 
alizing synthetic SO coupling in cold atoms 9 14^ opens 
intriguing new avenues to probe topological orders |15l - 
[26] with clean platforms in a fully controllable fashion. 
By far the experimentally realized SO interaction [10- 
[H] is the ID equal Rashba-Dresselhaus-type SO term 
through a two-photon Raman process as theoretically 
proposed in Ref. [9]. Detailed investigations show that 
only ID SO term can be realized in an atom gas with in- 
ternal A-type configuration as considered in experiments 
PI2Z1I1E]- With such ID SO interaction it is unfortu- 
nate not optimistic to reach the topological superfluid 
(SF) phase from an s-wave pairing state in cold atoms. 
The reason is because on one hand in a quasi-lD system 
no long-range s-wave order can be obtained; on the other 
hand, the proximity effect, as used in solid state exper- 



iments [T, is not realistic for cold atoms. Therefore, to 
find out a truly experimental scheme to observe MZMs 
with higher dimensional SO interactions is a foremost 
outstanding goal in the current field of cold atoms. 

In this letter, we propose a novel scheme to observe 
both the quantum anomalous Hall effect (QAHE) and 
chiral topological SF phase with (pseudo)spin-l/2 cold 
fermions trapped in a blue-detuned square optical lattice 
and coupled to two periodic Raman fields. The remark- 
able advantage of our model in the realization is that the 
square optical lattice and Raman fields, used to induce 
2D SO interactions, are generated simultaneously by the 
same laser fields. Basing our model on real cold atom 
systems and conventional technologies, we show that the 
predicted exotic topological orders can be studied with 
currently available experimental platforms. 

Model.— We start with the quasi-2D cold fermions 
trapped in a conventional square optical lattice, with 
their internal degree of freedom (atomic spins) experienc- 
ing periodic Raman fields M^^y induced by two-photon 
processes, as illustrated in Fig. [ij The dynamics of 
fermions are governed by the following Hamiltonian 
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where M^ = Mq sin(fcox) and Aly — Alosm^kox), and 
the optical dipole potential V{r) — —Vo[cos'^ (kox) + 
cos^(A;oy)] forms the square lattice. Below we describe 
how to obtain H with the realistic cold atom plat- 
form. The diagram for light-atom couplings is shown 
in Fig. lija). The ground states \g^,i) are coupled to 
the excited one |e) through two pairs of lasers rii^2 
and r2i^2, which form two A-type configurations with 
large one-photon blue detunings Ai and A2, respectively. 
For the present purpose, we set Oi = VIq sin(fcoa;) and 
^2 = ^osin(fcoj/), while fii = i(l2 — i^o to be constants. 
For the parameter regime that A12 ^ |Ai — A2I ^ 






1 


k' 


L J 


L 




^ ^-. 














^-- ja„^^ 


\. 


z^m 




r " 


r» T T J 


r 



-uy 



(b) 



FIG. 1: (Color online) (a) A realistic optical-dipole transition 
diagram in cold fermions (^''K) coupled to two pairs of laser 
beams r2i,2 and 17i,2 under large blue-detuning condition, (b) 
With this configuration a square lattice potential and two 
periodic Raman fields M^^y are simultaneously generated. 



|J7q/Ai_2|, the population of the excited state |e) is neg- 
ligible [S] and the hght-atom couphngs through two de- 
tuned A-type configurations lead to two independent Ra- 
man transitions between \g^^) and jg^), with the Raman 
fields given by M^ — h\^Q^ svaikQx) j ^\ and iMy — 
jfi|riopsin(/i;oy)/A2, respectively [Fig. [ija)]. Such two 
Raman processes can be precisely tuned with acoustic- 
optic modulator to have a small two-photon off-resonance 
(5 (|(5| <C |rio/Ai.2|) which gives rise to the Zeeman term 
?7i^ = hSl2 along z axis in H [T5] . Finally, note in the re- 
alistic atom system, e.g. for the ^°K atoms the lasers rii.2 
(f2i.2) also couple to the state \g\^ (l5t)) with correspond- 
ing large detunings Aqj ^ (Abia) [see Fig. IT] (a)]. While 
these additional couplings cannot lead to new Raman 
transitions between \g^,x) and other ground levels due 
to large two-photon detunings [HJ [13] , they contribute 
to the diagonal optical potentials for the two ground 
states and we get Vf = Ej=i,2(I^jP/^j + l^jlV^bJ 
and V^ = Ej(l%P/^j + miV^Qj)- In the realistic ex- 
periment (e.g. for '*'^K atoms) the magnitudes of detun- 
ings ^a,b,j (can be over 10"^ GHz) are much larger than 
their differences (in the order of lOMHz) [T^]. There- 
fore in the formulas of V^ and M^^y we can take that 
Ai = A2 = Aqj 2 = ^f)i,2 — ^; which is followed by 
Vf = Vi = -(|rJo|^/A)[cos^(fcoa;) + cos^(fco2/)] + const. 
and Mq = h\ilQ\'^/A. Neglecting the constant terms 
yields the effective Hamiltonian in Eq. (IT]). In brief, in 
Fig. [Tithe lasers fli,2 (red lines) generate the conventional 
square lattice. Together with another pair of lasers fli,2 
the Raman fields are simultaneously created. 

Before proceeding further, we provide two important 
remarks on the realization. First, since the Raman 
fields and the optical lattice potentials are simultane- 
ously achieved by the same lasers, the phase fluctuations 
in the lasers, characterized by fii = fig sin(fcoa; -I- ^i""^) 
and ^2 = ^0 sm{koy + ^f"'^), only lead to global shift of 
the lattice and Raman fields as illustrated in Fig. fljb) . 
The relative spatial profile of M^^y and V{r) is always 
automatically fixed and therefore the effective Hamil- 



tonian (IT]) is unchanged. This greatly simplifies the 
experimental set-up and makes the realization of this 
model is no more difficult than that of a conventional 
square lattice for spin-1/2 cold fermions. Second, the 
blue detuning is essential for the realization. If us- 
ing red-detuned lasers, i.e. ^a.bj < 0, one has then 
Vf^l = -(|rj^/A|)[sin^(A:oa::) + sin^(fco2/)], which shifts 
1/2 lattice site relative to the Raman fields. The red- 
detuned laser couplings cannot lead to 2D SO interaction 
or nontrivial topological orders as studied below. 

Quantum anomalous Hall effect.— The tight-binding 
model of H can be derived straightforwardly. We 
take that fermions occupy the lowest s-orbitals 4>scr 
{a =t, i), and consider only the nearest-neighbor 
hopping terms. The tight-binding Hamiltonian is 
given by i?Ti - -tsT,<lj>,a^t^ja + ^^^^(^t " 
%i)+[l]<?J>^lo4f%i +H-C-], where ts denotes spin- 
conserved hopping, the 2D lattice-site index i = (ix,iy), 
and n^ — & <> . From the even-parity of the s- 
orbitals and the profile of Raman fields shown in 
Fig. ITJb), it can be directly verified that the spin-flip 
hopping terms due to the Raman fields satisfy i^o'-'^^^ = 
±(-l)^'ii? and tlv'^^^ = ±i{-l)HfJ, where ti? = 
Mo J (Pr<f)s{x,y)sm{kox)4>s{x — a,y) with a the lattice 
constant. Redefining the spin-down operator c^, — )■ 
gi7rrj/a- ^ .^g rccast the Hamiltonian into 
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-[E ^- (lt^^==+u - 4t^^-u) + H-C-] + 

'[E*^-(4A+U"SV..-u)+H.c.]. (2) 



It is convenient to rewrite H in the fc-space and i^Ti = 



,(0) 



(0) . 



dx = 2tso sin(fcya), dy — 2tso sin(A:2;a) and d^ — —niz + 
2is cos{kxa) + 2tg cos{kya). A simple analysis shows that 
the bulk system is gapped when |mz| 7^ 4ts,0. By cal- 
culating the first Chern number C^^ — (47r)^^ J d^kn • 
3fc^n X dk n, with n = {dx,dy,dz)/\d{'k)\, one can ver- 
ify that the above Hamiltonian describes a topological 
insulator for < \mz\ < Atg and otherwise a trivial in- 
sulator. This result can be read by analyzing the band 
structure around the four independent Dirac points at 
kc = (0, 0), (0,7r), (7r,0), and (7r,7r), where the system 
are described by massive Dirac equations. First, it is 
straightforward to know that the system is a trivial in- 
sulator when m^ 3> 4is. Upon reducing rriz to m^ = 4is 
the band gap closes at kc — (0, 0) and reopens when 
< TO2 < 4ts , with the effective mass at this Dirac point 
changes sign and leads to the change by 1 in the Chern 
number. Therefore the system is in the topological phase 
when < ni^ < Atg with C^^ = 1. In the same way. 



one can determine that the phase is also topological for 
— 4is < niz < 0. On the other hand, both masses of 
the Dirac Hamiltonians around kc — (0, tt) and (tt, 0) 
change sign when varying m^ < to m^ > 0, implying 
that Cf'^irriz = 0+) - C^^im^ = 0_) = 2 and therefore 
C^^ = -1 for -4ts < niz < 0. Wc then conclude that 
the Hamiltonian (pi) realizes a QAHE with the quantized 
Hall conductance 



I'^AI 



E ft 



,AH 



sgn{mz)^, for < Im^l < 4is, 
0, for |m^| >4is,TO^ =0. 



(3) 



It is noteworthy that while the QAHE was previously 
studied in spinless cold fermion models j,l6 , which are 
very hard to be achieved in the realistic experiments, we 
have proposed here a novel scheme to observe QAHE in 
spin- 1/2 fermion system which is truly accessible with 
current experimental platforms. 

Chiral topological superfluid phase.— We next study the 
topological chiral superfluid phase by considering an s- 
wave interaction which can be well controlled with Fes- 
hbach resonance in cold atoms 29J. In optical lattice, 
this interaction is described by the attractive Fermi Hub- 
bard model Hint = ~ S?^"'?t'^4' '^'^^^^ the regular- 
ized 2D effective interaction U = h?asy/%TT/{mlzl\^) with 
h. = {h/riiuJzY''^ , hd — {h/'muj2dY''^ , ^'^^ the s-wave 
scattering length a^ > 0. The parameters W2d and lo^ 
denote the 2D square lattice trapping frequency and a 
tight trapping frequency along z axis, respectively. We 
first study the superfluid phase with the self-consistent 
mean fleld approach by introducing the s-wave superfluid 
order parameter by Ag = (C^/-/Vo) X]k;(^ktC-k^) with A^o 
the number of lattice sites. Then the total Hamiltonian 
in the Nambu basis ?Ak — (cktiCk^: cl^., , — clj^.^) can be 

written as -ffBdG = Z^k V'I'HBdG(k)V'k, where 



■HbcIG = dxCTx ®Tz+ dyGy ® Tz 

-f(Ag/(g)T+-|-H.c.). 
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Here ^ is the chemical potential, T^.y^z are Pauli matrices 
acting on the Nambu space, and t± — [t^ ± zTj,)/2. The 
s-wave order parameter is solved by the gap equation 
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FIG. 2: (Color online) Self-consistent calculations of the s- 
wave order As [for (a)] and the bulk gap Eg [for (b)] versus 
rriz and fi at zero temperature. Topological phase transition 
occurs at Eg = 0. The hopping coefficients 4o = is- 



topological insulator in proximity to an s-wave SC [1], 
which renders the 2D chiral p + ip topological SC phase. 
We therefore expect that in our model the chiral Ma- 
jorana edge modes and MZM in the vortex core can be 
obtained. Similar as the case in QAHE, the SF phase dia- 
gram can also be determined by analyzing the properties 
of the bulk gap. By examining E^ with rriz > we flnd 
that the bulk gap closes at k = (0, 0) and k = (tt, tt) for 
fi'^ + \As\'^ — {niz — 4ts)^ and (m^ -f 4ts)^, respectively. 
Furthermore, at the points k ~ (0, tt) and k = (tt, 0) 
the bulk gap closes if /x^ -I- jA^j^ — m?^. These prop- 
erties lead to a rich phase diagram as presented below. 
From the low-energy Hamiltonian T-Lbag we know that 
the superfluid phase with |As| > 0, ttt-z — 4is, and 
/i = is topological. This implies that if jm^j > 2ts, 
the SF phase is topologically nontrivial with Chern num- 



-1-1 for (JTO^I - AtsY < A*^ 



I A, 



< m 



where the energy spectra E^ = (F^ -I- jdkj^ ± 2wk)^/^ by Cf^ ^ -1 when (jm^ 



ber Cf^ 

and Cf^ = -1 for ml < /i^ + jA^I^ < (jm^j + Ats)'^. 
The change by 2 in the Chern number is originated from 
the fact that tuning /i^ -t- jAgj-^ from less than to larger 
than ml reverses the mass terms at both k = (0, tt) 
and k = (7r,0). On the other hand, when niz — the 
phase is always trivial. Actually, when tti^ = the low- 
energy Hamiltonian Hedo for the SF is captured by two 
Dirac cones around k = (0,7r) and k — (7r,0) with op- 
posite chiralities and having s-wave pairing. Such two 
Dirac cones resemble p + ip and p — ip superfluids, re- 
spectively, which couple to each other and cancel out, 
yielding a trivial phase. With this result we know that 
if \mz\ < 2ts, the SF phase is topologically trivial for 
ml < ^^ + |As|^ < (|m2|—4ts)^, while it is nontrivial with 

4t,)2. 



-4i.)2<^2_ 



< mz + 



diagonalizing HedG, with w^ == [fJ-'^idl + dl) + dlV^^/'^ 
and r2 = ^2 _^ |A,|2, and /? = l/iksT). 

To identify an effective p + ip chiral superfluid phase 
from the Hamiltonian "HedG: we analyze the special sit- 
uation with TTiz — -its- In this case the low energy spec- 
trum is captured by the Hamiltonian (up to linear or- 
der of the momentum) "HedG = 2tio {k^ay + kyax)Tz — 
^Tz + (AsT+ -l-H.c). It is interesting that this low-energy 
model is equivalent to the surface Hamiltonian of a 3D 



In cold atoms the parameters such as ttt,^, ts and // can be 
precisely adjusted in the experiment to tune the system 
into topological chiral SF phases. 

The self-consistent solutions for As and mean-field 
phase diagram are shown in Fig. [2] (a) and (b), respec- 
tively. It can be seen that the s-wave order As vanishes 
when the noninteracting Hamiltonian iJxi has a large in- 
sulating gap and the chemical potential /x locates deep 
in the band gap [left most area in (a)]. On the other 



hand, an appreciable As is obtained in the area with 
1.5fs < fj, < 3.5ts and < m^ < 3.0. In such param- 
eter regime the SO terms of dx and dy dominates over 
dz and a relatively large density of states at Fermi level 
is present. In Fig.l2](b) the superfluid phases with both 
Cf^ = +1 and —1 are obtained, respectively leading to 
chiral and anti-chiral Majorana edge modes localized in 
the boundary [T1. 

Note that in 2D superfluids no long-range order ex- 
ists at finite temperature and the critical temperature 
by mean-field theory is often overestimated. Instead, the 
Berezinsky-Kosterlitz-Thouless (BKT) transition occurs 
at the critical temperature which is limited by entropi- 
cally driven vortex and antivortex proliferation |30j . The 
BKT temperature is calculated by 



Tbkt — ^PsI^stTekt), 



(6) 



where ps is the superfluid stiffness (superfluid density). 
With the presence of phase fluctuation, the superfluid 
order parameter can be described by A^ = Aoe'^'-''', 
where 9 varies slowly in the position space. The sim- 
plest way to estimate the stiffness is to use its relation to 
the supercurrent density j^. — psV6'(r). For the case that 
9 varies slowly with position, we can approximate that 
q = V6'(r), which implies that the s-wave pairing occurs 
between two fermions with a center-of-mass momentum 
q. The current of the system can also be calculated by 
js = ^Tr[e-^^BdGai/BdG(M,Ao,V0)/aq]/Tr[e-'3^BdG]^ 

with A the area of the square lattice. This leads to 



Ps = 



I TT[e-P"^^^d^HBdG{l^,Ao,V9)/d{V9f 



qA 



Tr[e-'5^BdG] 



(7) 



Together with the gap equation for Ag we get the super- 
fluid stiffness. It can be verified that the above results are 
consistent with those through a standard functional path 
integral approach, where the phase fluctuation of the su- 
perfluid order parameter is described by an effective ac- 
tion Sfinc = {l/2)Jd^rp,i\79f. The BKT temperature 
is then given by Eq. (rol). We have numerically confirmed 
that when A^ = the supercurrent j^ vanishes. 

The numerical results for Tbkt versus U and p, are 
shown in Fig. [3] For the case that the Fermi energy 
touches the bulk edge, the BKT temperature increases 
gradually with U [blue line in Fig. IsFa)], while when the 
Fermi energy locates deep into the bulk, the critical tem- 
perature increases sharply with U and approaches a sat- 
urate magnitude in a small Hubbard interaction (red, 
green, and black curves). The BKT temperature versus 
II with fixed U = bts and ruz is plotted in Fig. ^h). It 
can be found that with rriz = Ms, the maximum Tbkt is 
around OMs- Note that the upper value of the blue de- 
tuning is limited by the hyperfine-structure splitting and 
can be taken as 2tt x f .7THz for ^°K atoms, which gives 
the recoil energy En/h ^ hk'^/2m = 27r x 8.5kHz using 
lasers of wavelength 764nm to form the square lattice. 
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FIG. 3: (Color online) BKT temperature as a function of 
\J with different values of m^ and /i (a), and versus p at 
\J — 5ts and with different magnitudes of rriz (b). The SO 
coupled hopping tL = is. Topological orders of different 
Chern numbers are indicated for the case with rUz — 3ts. 



Taking that f^o = 27r x 0.27GHz, we have Vq = 5Er, 
the lattice trapping frequency uj — 2TT x 34.9kHz, and 
the hopping coefficients tso /h ^ ts/h ~ 2tt x 0.52kHz. 
With this parameter regime we find that the bulk gap 
for QAHE Eg = 2tt x 1. 04kHz by setting ^ = 2*^, corre- 
sponding to the temperature T = 90nK for observation, 
and the maximum BKT temperature for the topological 
superfluid phase with U = 5is is about f OnK. 

Conclusions and Discussions.— We have proposed a 
truly experimental scheme to observe 2D SO interac- 
tion which leads to both the quantum anomalous Hall 
effect and chiral topological superfluid phase in spin-f /2 
cold fermions trapped with a blue-detuned square op- 
tical lattice and coupled to two periodic Raman fields. 
We calculated the bulk gap of the topological states and 
the BKT temperature of the topological superfluid phase, 
and show that the predicted topological orders are reach- 
able with the realistic parameter regimes. 

The present model exhibits many essential advantages 
in the realization. First of all, the square optical lat- 
tice and the Raman fields, used to induce 2D SO inter- 
actions, are generated simultaneously through the same 
laser fields, which greatly simplifies the setup for exper- 
imental studies. On the other hand, for blue-detuned 
optical lattice the atoms are trapped in the minimums of 
lattice potentials which may generally minimize heating 
effects |3T]. Note that while the topological superfluid 
phase with Majorana modes requires a fermion system, 
the topological insulating state can be demonstrated with 
cold bosons [35] , which makes experimental observations 
even more straightforward. It is worthwhile to point out 
that 2D chiral topological orders are classified by a Z in- 
variant, and is stable when stacking multi-layers of the 
2D system. This enables a study of topological superflu- 
ids in the multi-layer square lattices which may have a 
higher transition temperature. The present work paves 
the way to observe 2D topological orders in cold atoms, 
and will lead to the ultimate observation of non-Abelian 
Majorana fermions with realistic cold atom platforms. 
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discussions on the experimental feasibility of this work. 
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